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Abstract 

We compute the hadron mass spectrum, the quark masses and the meson decay 
constants in quenched lattice QCD with non-perturbatively 0{a) improved Wilson 
fermions. The calculations are done for two values of the coupling constant, /3 = 6.0 
and 6.2, and the results are compared with the predictions of ordinary Wilson fermions. 
We find that the improved action reduces lattice artifacts as expected. 
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1 Introduction 



The calculation of hadron masses in lattice gauge theory has a long history. Over the years 
there has been a steady improvement in computing power and methods allowing simulations 
on larger volumes and smaller quark masses with higher statistics. However, the progress to 
smaller lattice spacing a has been slower because of the high computer cost which increases 
by at least a factor of (1/a)^. Since it is so expensive to reduce cut-off effects by reducing a, 
we should consider reducing them by improving the action. 

A systematic improvement program reducing the cut-off errors order by order in a has 
been proposed by Symanzik [|I| and developed for on-shell quantities in ref. 0. The standard 
gluonic action has discretization errors of O(a^), but those for Wilson fermions are of 0(a). 
Therefore it is the fermionic action which is most in need of improvement. 

Sheikholeslami and Wohlert proposed the action (we assume r = 1 throughout the 
paper) 

where is the original Wilson action and 



F,Ax) = ^ E (Uix),, - Uix)l). (2) 



In eq. (0) the sum extends over the four plaquettes in the /iz/-plane which have x as one corner, 
and the plaquette operators U{x)^i, are the products of the four link matrices comprising the 
plaquettes taken in a clockwise sense. If csw is appropriately chosen, this action removes 
all 0{a) errors from on-shell quantities such as the hadron masses. A non-perturbative 
evaluation of this function leads to |Q 

1 - 0.656^2 _ 0.152^4 _ 0.054/ . 
cswig) = 1_ 0.922 ,2 ' ^ - 

When we talk about improved fermions in the following, we always understand that csw has 
been chosen according to eq. (^. 

In this paper we shall present results for the light hadron mass spectrum, the light and 
strange quark masses and the light meson decay constants using the improved action. The 
calculation is done for two values of the coupling, f3 = 6.0 and 6.2, which allows us to test 
for scaling. 

The mass calculations extend our earlier work where we have examined the csw 
dependence a.t (3 = 6.0. To exhibit the effect of improvement, we have also done calculations 
with Wilson fermions on the same lattices. Most of our Wilson data come from our structure 
function calculations p, , and we combine this with masses from the literature at other (3 
values to see the dependence on a clearly. 
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From the meson correlation functions we also extract meson decay constants and quark 
masses. However, simply improving the action is not sufficient to remove all 0{a) errors 
from these quantities. Here we also have to improve the operators which is done by adding 
higher dimensional terms with the same quantum numbers in an appropriate fashion. 

This paper is organized as follows. In sec. 2 we briefly describe our numerical method. 
The hadron masses are given in sec. 3, concentrating in particular on the extrapolation to 
the chiral limit and the scaling behavior of improved and Wilson action results. In sec. 4 
we compute the light and strange quark masses using two different methods, from the axial 
vector current Ward identity and from the lattice bare quark masses. The meson decay 
constants are discussed in sec. 5. Finally, in sec. 6 we give our conclusions. 



2 Computational Details 

Our calculations have mainly been done at (3 = 6.0 and 6.2 on 16^32, 24^32 and 24^48 
lattices. We use Quadrics (formerly called APE) parallel computers. For the improved case 
the parameter csw is given from eq. (|D as csw = 1-769 at /5 = 6.0 and csw = 1-614 at 
/3 = 6.2. The simulations are done for at least five different k values in each case. This helps 
with the extrapolation to the chiral limit. 

For the gauge field update we use a combination of 16 overrelaxation sweeps followed 
by a three-hit Metropolis update. This procedure is repeated 50 times to generate a new 
configuration. 

The improvement term in eq. (|T|) appears in the site-diagonal part of the action. The 
major overhead in our case is multiplication by this term during inversion of the fermion 
mass matrix. In our basis of hermitean gamma matrices we can rewrite this term as P| 

1 - ^K,gcswo- ■ F = 




1 1 
-1 1 



(4) 



where A, B are 6x6 matrices (two-spinors with color), so that instead of a 12 x 12 multi- 
plication we have two 6x6 multiplications and two inexpensive coordinate transformations. 
This reduces the overhead for the improvement in the inverter from 45% to 30%. Also, the 
inverse of the matrix in eq. is required on half the lattice due to the even-odd precondi- 
tioning. We now have to invert two 6x6 instead of a 12 x 12 matrix. However, this is only 
required once for each propagator inversion. 

For the matrix inversion we mainly used the minimal residue algorithm, except for the 
lightest quark mass on the larger lattices where we used the BiCGstab algorithm |TU|. As 
convergence criterion we chose 

|r| < 10"^ (5) 
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for the residue, which is the best that can be achieved for our single precision machine. 
For the mass calculations we used Jacobi smearing for source and sink. For a detailed 



description of our application of this procedure see ref. |TT[. We have two parameters we 
can use to set the size of our source, the number of smearing steps, Ng, and the smearing 
hopping parameter, k^. We chose Ng = 50 for (3 = 6.0 and 100 for (3 = 6.2 and Kg = 0.21 
at both (3 values. This gives roughly the same r.m.s. radius in physical units in both cases, 
namely 0.4 fm. To define the matrix elements for the decay constants and quark masses, we 
have also computed correlation functions with smeared source and local sink. This does not 
require any additional matrix inversions. 

At P = 6.0 and csw = we had generated 0(5000) configurations for our structure 
function project on which we have computed the hadron masses. To these we added 0(150) 
new configurations on which we computed the meson decay constants and the chiral Ward 
identity. For csw = 1-769 we have analyzed 0(1000) configurations. For the heavier quark 
masses, k, = 0.1487 and k = 0.1300, 0.1310, 0.1320, respectively, the number of configurations 
was 0(200). On the 24'^ lattice we have generated O(IOO) and 0(200) configurations at 
Csw = and 1.769, respectively. At (3 = 6.2 we only ran on 24^ lattices. Here we have 
analyzed O(IOO) configurations for csw = and 0(300) configurations for csw = 1-614. 
We employed both relativistic and non-relativistic wave functions |^, except for the high 
statistics runs where we only looked at the non-relativistic wave function in order to save 
computer time. 

Besides our calculations at /5 = 6.0 and 6.2 we also made exploratory studies at /5 = 5.7 
to see what effect varying csw has on coarser lattices. If one decreases P, increases csw or 
increases k, one starts to get problems with exceptional configurations. This showed up in 
non-convergence of our fermion matrix inversions. It was, however, only a real problem at 
(3 = 5.7, Csw = 2.25 and § /3 = 6.0, csw = 3.0. 



3 Hadron Masses 

We consider hadrons where all the quarks have degenerate masses. We looked at vr, p, 
nucleon (A^), oq, ai and bi masses, and we have used this nomenclature for all quark masses, 
not just in the chiral limit. 

In our mass calculations we have made single exponential fits to meson and baryon cor- 
relators over appropriate fit ranges. The errors are determined using the bootstrap method 
with 50 data samples. We present our hadron mass results in tables ^ |] and ^ Table ^ 
updates the results presented in ref. For the meson masses we found very little difference 
between using relativistic and non-relativistic wave functions, and we settled for relativis- 
tic wave functions (except for the high statistics runs). For the nucleon we have chosen 
non-relativistic wave functions which performed slightly better because the effective mass 
plateaus extended to larger times. At (3 = 6.0 we repeated the lightest quark mass on 16^32 
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on the 24^32 lattice, for both improved and Wilson fermions. The values agree within less 
than 3%. This indicates that all our results on the 16^32 lattice do not suffer from significant 
finite size effects. 



Chiral Behavior 

To obtain the critical value of k, Kc, and the hadron masses in the chiral limit, we extrapolate 
our data to zero vr mass. We first tried 

Using this relation gives a rather poor fit of the data, and we saw that there was a slight 
curvature in a plot of against Quenched chiral perturbation theory predicts 



m 



2 



1 1 \ — 



where 5 is small and positive. We made fits using this formula but found that 5 was always 
negative. As in our previous work we conclude that our k values are too far from for the 



formula to be applicable. This is in agreement with observations made by other authors |T3 
As an alternative parameterization of the curvature we used the phenomenological fit 

1 1 

- = — + h2ml + h^ml. (8) 

In table ^ we give the values of for the different fits. The linear fits give x^/dof values of 
up to 40. The other two fits both give acceptable values of x^, but eq. (||) usually gives a 
lower than eq. (|^. In the following we shall take Kc from the phenomenological fits. 

In fig. Ill we plot for improved Wilson fermions. We compare our results with the 
results of ref. The agreement is excellent. In one-loop perturbation theory is given 
by i 

= i[l + /(0.108571 - 0.028989C5W/ - 0.012064 4^^)]. (9) 
8 

The tadpole improved value of Kc that follows from this result is 

= 1[1 + ^2(^0.025238 - 0.028989 C5VI/M0 - 0.012064 {cswulf)]uo\ (10) 
8 

where csw is given by eq. (0), 

no = (^Trf/o)3 (11) 
and g*^ is the boosted coupling constant defined by 

9*' =9^/4. (12) 
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In fig. ID we compare the tadpole improved perturbative formula (|T^) with the data where 
for the larger couplings we have taken uq from [H, T^]. The curve and the data points 
agree within less than 1%. In eq. (|10]) one has the choice of using the lowest order tadpole 
improved value of csw, namely Uq^ [|[, or the value from eq. which is the value actually 
used in the simulations. Both procedures remove all the tadpole diagrams and differ only 
by small 0{g^) terms, so they are both reasonable. We prefer the second choice. 

We fit the other hadron masses by the formula 

mjj = bo + hml + h^ml , H = p,N,---. (13) 

The result of the fit is shown in fig. ^ for both improved and Wilson fermion data. The 
Wilson fermion data are the world data compiled in tables ^, || and 0. 

We find this to be a more appropriate fit formula than the ansatz I^Gj 

ruH = fe'o + b'^ml + b'^ml, (14) 
because for the nucleon the plot of against (or 1/k) is less curved than rriN against 



m: 
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(Note that the two formulae differ only by terms of 0(m^)). To decide which fit formula 
is best and to do a reliable extrapolation to the chiral limit, it is important to have many 
K values. For the oq, ai and bi masses only a two-parameter fit with 63 set to zero was 
reasonable. The mass values in the chiral limit for our data are also given in tables |1|, |^ and 

I- 

We see that the effect of improvement is largest for the p mass. In the chiral limit the 
difference between improved and Wilson results is 25% at (3 = 6.0 and still 12% at [3 = 6.2. 
It is quite common to define the physical scale from the p mass. The relatively large change 
of this quantity from the Wilson to the improved case suggests that it contains large 0(a) 
corrections, and that this procedure is misleading. A better procedure is to use the string 
tension or tq the force parameter, as the scale. For the nucleon mass the difference 
between the two actions is smaller. 



APE Plots 

In figs. ^ and ^ we show the dimensionless ratio 771^/ nip as a function of {mT^/mpY , a so- 
called APE plot, for (3 = 6.0 and 6.2, both for improved and Wilson fermions (the latter 
using the world data given in tables ^ and |^). The solid lines are the results of the ratio of 
the fits in fig. ^ At /? = 6.0 we find that the mass ratio data are rather different for the two 
actions. The improved results lie consistently lower than the Wilson results. At (3 = 6.2 we 
find the same pattern in the data. 

At (3 = 6.0 we can say something about the chiral limit. Our fits give m^/rrip = 1.20(6) 
for improved fermions and mj^/mp = 1.33(2) for Wilson fermions. The improved results 
come closer to the physical value than the Wilson results. At (3 = 6.2 we are lacking data at 
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small quark masses and on larger volumes. In the chiral limit our fits give m]si/mp = 1.32(11) 
for improved fermions and mj^/mp = 1.39(12) for Wilson fermions, so that we cannot say 
anything conclusive about the behavior of the two actions in the chiral limit in this case. 



Scaling Behavior 

Let us now look at and compare the scaling behavior of the two actions. We shall limit 
our discussion to the p mass because the errors of the nucleon are too large to make precise 



statements. In order to exhibit the cut-off effects most clearly, it has been suggested [|28[ that 
rrip should be plotted in units of the square-root of the string tension K which has cut-off 
errors of O(a^) only. In table |]we have compiled the world string tension data. When there 
are several calculations, we performed the weighted average. 

In fig. ^ we plot the ratio mp/ \/K as a function of a\/K. This is done for fixed physical 
TT masses with = 0, 2K and 4K. Comparing hadron masses at larger quark masses has 
the advantage that this does not require large extrapolations of the lattice data but rather 
involves small interpolations only. The Wilson fermion data shown are a fit to the world 
data compiled in tables |^, ^ and |^. As expected, the Wilson masses show practically a linear 
behavior in the lattice spacing a. We have done a simultaneous linear plus quadratic fit to the 
Wilson data and a quadratic fit to the improved data. The fit is constrained to agree in the 
continuum limit. The result of the fit is shown by the solid lines in fig. ^ In the continuum 
limit we obtain nip/y/K = 1.80(10). We compare this result with the experimental p mass. 
For the string tension we take the value 



K = A27MeV (15) 



which has been obtained from a potential fit to the charmonium mass spectrum |^6[. Using 
this value the physical nip/^/K is 1.80 which agrees with the lattice number. 

As mentioned previously, an alternative scale from the potential is tq. We have also 
compiled lattice results for tq in table |. We see that it scales very well with \/K, as the 
product tq^/K is approximately constant at about 1.19, while the lattice spacing a changes 
by a factor of more than five. However, the physical value of rgy/K is 1.06, taking r^^ as 
402 MeV which follows from the same potential that gives \/K = 427 MeV |^ . It does not 
seem that this discrepancy will vanish as a — > 0. It is telling us that the lattice potential has 



a slightly different shape to the continuum potential. This may be an effect of quenching |37 

Although at /5 = 5.7 we do not know the correct value of csw: using our larger value 
Csw = 2.25 we find nip/y/K = 1.94 in the chiral limit. Comparing this number with fig. 
it indicates that O(a^) effects are moderate even at this coupling. 
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Mass Splitting 



The vector-pseudoscalar mass splitting 

Av-ps = my - rrips (16) 

is experimentally rather constant for all quark flavors. One finds 

ml -ml = 0.57 GeV^ 
m^. -m^ = 0.55 GeV^ (17) 
-ml = 0.55 GeVl 

Quenched lattice calculations with Wilson fermions are unable to reproduce these numbers. 
Wilson fermions give a splitting which is much too small. In fig. |^ we compare the experi- 
mental values of — and mj^, — m\ with the lattice data and the mass fits. As before, 
we have taken the string tension eq. (|l^) as the scale. In fig. ^ we also show the results for 
improved fermions and the corresponding mass fits as well. There is a noticeable change 
when going to the improved case. We find good agreement with experiment for the absolute 
values. 

In the heavy quark effective theory |]38 



Ay_P5(x (f'a^.F^,*), (18) 

where \1/ is the heavy quark field. So it is natural that turning on the Sheikholeslami-Wohlert 
term would increase the mass splitting, and this is what we see. 



Wilson Kc 

Let us now come back to the critical value of k for Wilson fermions. In table ^ we have given 
the values of from a fit of the world data in tables ||, |^ and using the phenomenological 
ansatz (H). In fig. ^ we plot these results as a function of aVK (the string tension being 
taken from table |^). We see that Kc is a linear function of a over the whole range of the 
data which extends from (3 = 5.7 to 6.4. Comparing this with the improved Hc, which is 
approximately constant, we conclude that the Wilson Kc has large 0{a) effects. We also 
compare the Wilson data with the predictions of tadpole improved perturbation theory as 
given by eq. ( p!0D with csw = 0. Here we have taken the one-loop perturbative formula for a 
beyond (3 = 6.8 where there are no numerical values for the string tension available any more. 
Not even at the smallest value of a can perturbation theory describe the Wilson data. For 
improved fermions, on the other hand, the agreement with tadpole improved perturbation 
theory is quite good, as we have already noticed. 
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4 Quark Masses 



We shall now turn to the calculation of the quark masses. When chiral symmetry is dynam- 
ically broken, care has to be taken in defining renormalized masses. In the continuum the 



renormalized quark mass at scale p = /i can be written 



^Tr[5/(m)-^^i(0)]=m(/i), (19) 

where Sp is the renormalized quark propagator which is to be evaluated in a given gauge. 
This definition refers to the momentum subtraction scheme. It is usual to give the quark 
masses in the MS scheme. To convert from one scheme to the other, one has to go to high 
enough scales so that one can use perturbation theory. If the quark mass is defined in this 
way, then the renormalized mass is proportional to the bare mass. 

On the lattice the standard assignment of the bare mass is 

am{a) = \{--^), (20) 

giving the renormalized mass as 

m^(/i) = Z'^iafi, am) m(a), (21) 

where Z^^{afi, am) is the mass renormalization constant. We call this method of determin- 
ing the renormalized mass the standard method. 

An alternative way of defining a bare mass is by means of the PC AC relation between the 
divergence of the axial vector current = ip'j^'y^ip the pseudoscalar density P = ip'y^ip, 

. d,{A^{x)0) 
"^^^^ = 2{Pix)0) ' 

where O is a suitable operator having zero three-momentum and no physical overlap with 
A^lx) and P{x) to avoid contact terms. (See later on for a precise definition.) All operators 
are bare operators. To avoid anomaly terms in eq. (^), fiavor non-singlet operators are 
taken. We call this method the Ward identity method. The renormalized mass is then given 
by 

m^ifi) = -Jd^^^m{a), (23) 
Zp^{afi, am) 

where Zyi(am) and Zp^{an, am) are the renormalization constants of the axial vector current 
and the pseudoscalar density, respectively. 

The quark mass inherits its scale dependence from the renormalization constants Z^ and 
Zp which involve logarithms of fi. In the following we will compute Zm and Zp perturbatively 
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to one-loop order for lack of a better, non-perturbative determination. To keep the logarithms 
under control it is best to take a/x = 1 and do the transformation to any other scale by the 
renormalization group formula 

_8_ 

In the continuum limit both procedures should give identical results for m*^'^(/i). Note, 
however, that the two bare masses m and m can be different, though they both vanish in 
the chiral limit. On the lattice the two procedures may give different results for m*^'^(/i) due 
to non-universal discretization errors. 

The lattice calculation of the quark masses now proceeds in two steps. In the first 
step one has to find the k values corresponding to the real world by adjusting (e.g.) the 
pseudoscalar meson masses to their experimental numbers. In case of the Ward identity 
method one furthermore has to compute rh. In the second step the bare quark masses have 
to be converted to renormalized masses. We shall compute the masses of the u and d quarks, 
which we assume to be equal, and the mass of the strange (s) quark. 



Improved Fermions 

Let us consider the case of improved fermions first. Later on we shall compare our results 
with the predictions of Wilson fermions to see the effect of improvement. 

We will discuss the Ward identity method first. For the operator O we take the pseu- 
doscalar density 

P{{]) = Y,P{x^ = Q,x) (25) 

X 

and smear it as we did in the hadron mass calculations. As the -P(O) part is common to 
all two-point functions, we could have used any operator projecting onto the pseudoscalar 
state. Similarly, we write 

A^{t) = Y,A^{xi = t,x). (26) 



For improved fermions the axial vector current in eq. ( P^D is to be replaced by 

^4 ^ ^4 + CAadiP{x), (27) 

where ca is a function of the coupling only. The time derivative is taken to be the average 
of the forward and backward derivative. The coefficient ca has been computed in giving 
Ca = —0.083 at /3 = 6.0 and ca = —0.037 aX (3 = 6.2. The resulting bare mass 

^^"^ = WWW) ^ ^ 
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has been plotted in fig. | for /3 = 6.0 and our smallest quark mass on tlie 16^32 lattice. 
In fig. ^ we show the same quantity for /? = 6.2 and our smallest quark mass on the 24^48 
lattice. (Also shown in these figures are the results for Wilson fermions which we will discuss 
later on.) Equation (^) should be independent of t, except where the operators physically 
overlap with the source, if the cut-off effects have been successfully removed. In both cases, 
but in particular at (3 = 6.2, we see a smaller deviation from the plateau at small and large 
t values. To obtain the mass, we fit the ratio ( p8D to a constant. We have used the same 
fit ranges as for the pion mass. The results of the fit are given in tables 0. At f3 = 6.0 
in the improved case we see that at k = 0.1342 we have small finite size effects, indicating 
again that our results on the 16^32 lattice are not significantly volume dependent. 

For both the Ward identity and the standard method we choose to determine the k 
values from the pseudoscalar meson masses. Sometimes the 0(1020) meson is taken for the 
determination of the strange quark mass. However, we do not think that this is a good idea 
because of potential uo — cf) mixing We generalize eq. (§) to the case of two different 
quark masses by writing 

\{— + —)-- = hmls + hmls (29) 

/ Kl K2 Kc 

with the same coefficients 62, &3 as before. This is inspired by chiral perturbation theory 
where it is expected that the pseudoscalar mass is a function of the sum of quark and 
antiquark mass, + mq, even when quark and antiquark have different flavors. By fixing 
mps to the physical pion mass m^ri, using the string tension values compiled in table ^ 
with eq. (|T5|) as the scale, we find the value for Ku^d = ki = k,2- The strange quark mass is 
obtained by identifying mps with the kaon mass mx±, taking ni = Hu,d as input and solving 
for K,2 = i^s- This gives for the light mass 

1 M _ ^\ _ r 0.001836(36) for P = 6.0, 
^•'^ 2 kJ \ 0.001384(36) for /3 = 6.2. ^""^^ 

For the strange mass we get 

ma--(-~-]-l 0-0419(11) for P = 6.0, 
^ 2 kJ \ 0.0310(11) for/? = 6.2, ^'''> 

where rriu^d = 1/2 (m„ + m^). 

The bare masses rhu,d, fhs are computed analogously. We write 

m = ^(rfii + 7712) = 62"T'P5 + &3"^P5- (32) 

Using this parameterization we first fit the masses in tables |TU|, [TT] to the pseudoscalar masses 
in tables ||, ^. This gives us 62, We then determine fhu^, rhs by fixing mps to the physical 
pion and kaon masses, respectively, as before. 
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The mass dependence of the renormahzation constant Z^(am) can be parameterized 



as 



Z^(am) = (1 + 6^am)Z^. 



(33) 



The renormahzation constant has been computed non-perturbatively in ref. [Q. The 
fit formula in this paper gives = 0.7924 at /5 = 6.0 and Za = 0.8089 at (3 = 6.2. The 
coefficient 6a is only known perturbatively to one- loop order [|^. The best we can do at 
present is to take the tadpole improved value. For the boosted coupling we use af^^{l/a), 
giving 



bA 



a 



MS/ 



1/a) 1.912, 



(34) 



where we take af ^(1/a) = 0.1981 at f3 = 6.0 and af ^(1/a) = 0.1774 at (3 = 6.2 fg. For 
Zp^{afi, am) we write 

Zf^{afi, am) = (1 + bpam) zf^ (aii) . (35) 
The renormahzation constant Zp^{afi) has been computed perturbatively |0. The result 



is 



Z^'^iafi) = 1 - ^— C^(-61n(a/i) + 22.595 - 2.249c5vk + 2.0364^^), 



(36) 



with Cp = 4/3. We shall take the scale /i = 1/a and use the tadpole improved value of 
eq. ( p6|) which turns out to be 



Z*^^(a/i = 1) 



4:71 



16.967 - 2.999cswUo + 2.715{cswUo 



3\2 



Uq. 



(37) 



(We use Mo = 0.8778 at /3 = 6.0 and uq = 0.8851 at /5 = 6.2). The coefficient bp has also been 



computed perturbatively to one- loop order ||4^. Again we shall use the tadpole improved 
value 



bp = l + af ^(1/a) 1.924. 



(3^ 



We have also computed the renormahzation constants ^^(am), Zp{aiJ,,am) non-pertur- 
batively [^. So far we have results for /3 = 6.0 only. Our numbers are in fair agreement with 
the non-perturbative calculation in ref. |^ and the tadpole improved value (|37|). However, 
for small /i the constant Zp behaves very differently from the perturbative formula. 

To compare the results at the two different (3 values, we rescale them both to /i' = 2 GeV 
using formula (0). As before, we use the string tension to convert the lattice spacing into 
physical units. The resulting quark masses m^^^{2 GeV), m^'^{2 GeV) are given in table |1^. 

Let us now discuss the standard method. We already have determined ,^(a), ms{a) in 
eqs. (^), (0). For the renormahzation constant Z^^(afj.,am) we write 



Z^^{a^, am) 



[l + bmam)Z^'^{an). 



(39) 
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The constant Z^^(an) has been computed perturbatively [0. We obtain 



Z^^a/j) = 1 - -j—CF{6\n{afi) - 12.952 - 7.738csw + l.S804w)- 
The tadpole improved value, which we will be using, is 



z^'iafi = i; 



a 



MS, 



1/a) 



4:71 



(-4.110 - 10.3l7cswul + lMO{cswul) 



(40) 



(41) 



The coefficient bm has been computed in [43 1. The tadpole improved value is 



.^-af^(l/a) 1.210. 



(42) 



Again we extrapolate the quark masses to fi' = 2 GeV using eq. (^4]) . The results which 



follow from this approach are listed in table |T2| as well. 

The results of the Ward identity and the standard method may differ by O(a^) effects, 
and they do. We can 'fit' the a dependence by 



"^f = Co + caa^ 



(43) 



The result of the fit is shown in figs. |TU], |TT]. The continuum values from this fit are given 



in table |T2|. We find that the two methods give consistent results in the continuum limit. 
Taking the statistical average of the two results we obtain the continuum values 



^^(2GeV) 
(2GeV) 



m 



MS 



5.1 ±0.2MeV, 
112±5MeV. 



(44) 
(45) 



The Ward identity method appears to have larger O(a^) effects than the standard method. 

We may compare our results with the prediction of chiral perturbation theory, which 
cannot give absolute values but can determine the ratio of rus to mu,d- A recent calculation 



gives [46| mjmu d = 24.4 ± 1.5. We find m,/m„ d = 22.2 ± 1.2. 



Wilson Fermions 



Let us now consider the case of Wilson fermions. We proceed in the same way as before. 
The situation here is that Z^(am) is known non-perturbatively only for j3 = 6.0 and 
that Ba, bp and bm are only known to tree-level. So for Za we use the tadpole improved 
perturbative value 



Za 



An 



7.901 



Mo, 



(46) 
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and for 6^, bp and we take the tree-level results. Comparing with the non-perturbative 
determination at = 6.0 as well as with a non-perturbative calculation at /3 = 5.9, 6.1 
and 6.3 using the Ward identity |^8|, we find good agreement. The renormalization constants 
(a/i = 1) and Z^^{afi = 1) are obtained from eqs. (|37D, Pl|) by setting csw = 0. The 



resulting quark masses are given in table 0, and they are plotted and compared with the 
improved results in figs. 0, [Tl|. In this case we expect discretization errors of 0(a) instead 



of O(a^). So it is not surprising that the Ward identity and the standard method give results 
which are far apart. We find that the Ward identity method gives mass values which are 
closer to the continuum result. 

Finally, in figs. |T2|, we compare our improved quark masses with the world data of 
Wilson quark masses as compiled in ref. |^ for the standard method. These authors use the 
p mass extrapolated to the chiral limit to set the scale. At (3 = 6.0 the scale set by the string 
tension and by the Wilson action p mass differ by about 20% which explains the difference 
between our Wilson data and the world data in figs, [l! 
improves the scaling behavior. 



13i We see that the improved action 



5 Decay Constants 

The pion decay constant is well known experimentally and can be determined from the 
two-point correlation functions on the lattice as well, allowing for a further test of scaling of 
the improved theory. We shall also look at the decay constants of the K, p, K* and the ai 
meson. 

In Euclidean space at zero three-momentum we define 

(0|A|vr) = m^/^, 
(0|A|ai,A) = e{X),mlj,,, (47) 

2 

Tfl 

(0|V,|p,A) = e(A),-^, 
Jp 

where A and V are the renormalized axial vector and vector current, respectively, and e(A) 
is the polarization vector with Z^a (A)ej(A) = 5ij. The pseudoscalar and vector states are 
normalized by 

{p\p') = {2t, f2p,5{p-p'). (48) 

Note that our fa^ is defined to be dimensionless. 

In the improved theory the renormalized operators are 

A = {l + hAam)ZA{A^ + CAad^P), (49) 
= {l + bvam)Zv{y^ + icvadxT^x), (50) 
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where = ip'y^ip and T^,^ = ipa^^ip are the vector and tensor operators, respectively. We 
use the definition o"^,^ = i[7^,7j,]/2. Both currents are (partially) conserved, and hence no 
scale enters into their definition. The renormalization constant and the improvement 
coefficients ca and Ba have already been given in the last section. The renormalization 



constant Zy and the coefficients by and cy have been computed non-perturbatively in ref. [42 



m. At (3 = 6.0 the values are Zy = 0.7780, by = 1.472 and cy = -0.32(6), and at /3 = 6.2 



the numbers are Zy = 0.7927, by = 1.409 and cy = —0.22(7). While for most of these 
quantities the authors have given fit formulae in g"^, for cy we have read the numbers from 



the graph in |^0[, as no such formula exists yet. We have also determined Zy and by at 
(3 = 6.0 from our nucleon three-point functions and find consistent results. 

On the lattice we extract the meson decay constant from two-point correlation functions. 
For large times we expect that 

Co^oAt) = {Oiit)Olm 

^ '{0\Oi\H){0\O2\H)*e-"'»' + {Q\O\\H)*{0\Ol\H)e-"'"^^-'^] (51) 



2mH 



where 0{t) is of the form Vs ^ Y^x^i^i i)i being the spatial volume of the lattice, 

and (9^ = T]0 with 77 = ±1 being given by 74r'''74 = rjT. The rj factor tells us how O behaves 
under time reversal, i.e. whether the two-point function is symmetric or antisymmetric with 
respect to t T — t. Here T is the temporal extent of the lattice. In general we have 
computed correlation functions with local (L) and smeared (S) operators. 

We shall now consider the appropriate matrix elements separately. We start with those 
matrix elements necessary for the vr. With our conventions we set 

(0|A4|7r) = m./f , 

(52) 

{Olad^Plir) = — sinh am7r(0|P|7r) = m^a/^-*, 

where f^^\ f^^^ are defined to be real and positive. By computing C^fp and Cpf, we find for 
the matrix element of ^4 from eq. ( ^2]) 



r(0) o V^^ttAj^ p 

m^fn = -2/t , (53) 

' ASS 

/ipp 

and for the matrix element of d^P we obtain from the ratio of the Cpf, and C^fp correlation 
functions 

^ = sinham^-^. (54) 



15 



Alternatively, we can take the time derivative from the plateau in the correlation function. 
Numerically we found that it made very little difference to the result. 

For the ai we set 



(0|A|ai,A) = e(A),m2jW, 



and we find 



For the p we set 



km: 



(0) 



2k 



(0|a94rj4|p, A) 



3 J2k ^if 

e(A).mJ/W, 



(55) 
(56) 



sinh amp (0 1 Tj4 1 p, A) = ie(A)j mpfl/'-"'^^ 



and we obtain 



and 



2k 



^2rap J2k ^VkVk 
3 Efc ^VkVk 



(0) 



— i sinh am, 



(57) 



(58) 



(59) 



In tables and |14| we give the lattice results for the matrix elements calculated from 
the above formulas. The fits to the correlation functions, as for the masses, are all made 
using the bootstrap method. 

Collecting all the terms, the physical decay constants are given by 



fn 
fai 



;i + 6^am)Z^(/f +CAa/«; 
'l + bAam)ZAfi% 



(60) 



When the improvement terms are weighted with the appropriate c factors, they contribute 
about 10-20% at /3 = 6.0 and up to 10% at /5 = 6.2. It is thus important to improve the 
operators as well. 

To perform the chiral extrapolation, we make fits similar to those for the hadron masses, 
namely 



J n 

f 

J ai 



bo + b2ml + b-sml, 
bo + b2ml + fogm^, 
bo + b2ml + b^ml. 



(61) 
(62) 
(63) 
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We decided to fit the square of tlie decay constants ratlier tlian tlie decay constants tliem- 
selves because tliis sliows less curvature. Tlie fits and the data are shown in fig. n for f^^ 
and fp. We compare this result with the meson decay constants computed with the Wilson 
action. These follow from eq. ( pO| ) with ca,cv = 0. For Za we use the tadpole improved 
value given in eq. (|^), and for Ba we take the tree-level result {Ba = !)• The renormalization 
constant Zy (in the chiral limit) has been determined non-perturbatively from a two-point 



correlation function of the local vector current |H8| at /3 = 5.9, 6.1 and 6.3. Unlike the case of 



Za, we find significant differences between this determination and our determination using 
the nucleon three-point function. The latter gives Zy = 0.651(15) aX j3 = 6.0 which is close 
to the tadpole improved result. This indicates large 0(a) effects. Since we are applying Zy 
to a two-point function, we chose to use the non-perturbative result from ref. We inter- 
polate this result to /5 = 6.0 and 6.2 and find Zy = 0.565 and 0.618, respectively. For by we 
again take the tree- level result. Although the individual contributions of the improvement 



terms are significant, the overall result for /^r in fig. |TJ is not much changed when com- 
pared with the Wilson case for smaller quark masses. For larger quark masses, especially at 
f3 = 6.0, the Wilson /^r is larger. The situation is different for fp. Here we find a systematic 
difference of 10-20% at /3 = 6.0 and approximately 10% at /5 = 6.2 for all quark masses. In 
both cases the difference between the two actions becomes smaller with increasing /3 as one 
would expect. 



Our results extrapolated to the chiral limit are given in table |T5|, and we compare /^r 
and fp with experiment in fig. 15. For f^^ we find reasonable agreement of the improved 
results with the experimental value using, as before, the string tension as the scale. When 



including the data of ref. one sees that the Wilson results he lower, and it appears that 
the values are increasing as we approach the continuum limit. For fp both our improved 
and Wilson results lie within 5% of the experimental value. There is, however, a definite 
difference as we previously remarked. The Wilson numbers lie above the experimental value, 
while the improved ones lie below. One must remember though that in the Wilson case 
there is a systematic error in the renormalization constant Zy which may be larger than 
the statistical errors in the figure. The experimental number for the decay constant of the 
Oi is 1^ = 0.17(2) (in our notation). The agreement between experimental and lattice 
values is encouraging. 

We can avoid errors from extrapolating to the chiral limit by considering quark masses 
within our data range, as we have already done in figs. |^ and ^. The most physical k, 
values to use are those corresponding to the K mass. To obtain the decay constants we 
take eqs. (0) and ( |63D at = rriK- (Remember that we are using as a generic name 
for the pseudoscalar meson mass.) We give the results for fx and fx* in table |15|, and in 
fig. 16 we show the scaling behavior together with the experimental value for fx- We find 
the errors to be substantially reduced. For fx we see no difference between improved and 
Wilson results, both lying 10% below the experimental value. For fx* the error bars have 
become small enough to attempt an extrapolation to the continuum limit. The curves are a 
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simultaneous fit, linear for the Wilson and quadratic for the improved data, constrained to 
agree in the continuum limit. In this quantity there appear to be large O(a^) effects in the 
improved case. 

6 Conclusions 

The goal of this paper was to investigate the scaling behavior of 0(a) improved fermions. 
If scaling is good, the results we get should already be close to the continuum values for 
present values of the coupling. To this end we have done simulations for two values of (3 and 
looked at two-point correlation functions from which we derive hadron masses, quark masses 
and meson decay constants. 

First we looked at hadron masses. The most visible difference between Wilson and 
improved fermions is that the p mass is much lighter in the Wilson case at comparable 
pion masses. In fig. ^ we see that the improved action has brought the p mass closer to 
its physical value when we use the string tension to set the scale. In this figure we have 
compared the Wilson action p masses at many different scales. We see a linear behavior in 
the lattice spacing a as one would expect. For improved fermions we find the discretization 
errors reduced for our couplings. 

A problem with Wilson fermions was that they could not describe the vector-pseudoscalar 
mass splitting adequately. This problem seems to be cured by using improved fermions. 

Quark masses are important parameters in the Standard Model. Experimentally, their 
values are poorly known, and a reliable lattice determination would be useful. Using two 
different methods, we have determined the light and strange quark masses. Our results can 
be seen in figs. [1^, 0. Both methods give consistent results for improved fermions. In the 
continuum limit we find for the average of u and d quark masses m^^^{2 GeV) = 5.1±0.2 MeV 

and mf^^{2GeV) = 112 ± 5 MeV. In the Wilson case the discrepancy between the two 
methods is much larger, hinting at substantial 0{a) effects. 

When calculating the decay constants, an advantage of using the improved theory is 
that the renormalization constants and improvement coefficients for /^j, fp and fx* are 
known. For fx we still have to use the perturbative values of La because they have not yet 
been computed non-perturbatively. A systematic uncertainty in the Wilson case lies in the 
choice of the renormalization constants. While the results are in reasonable agreement with 
phenomenology, the data are at present not precise enough to discuss an extrapolation to 
the continuum limit, with the possible exception of Jk*. In that case it looks that there are 
relatively large O(a^) effects between (3 = 6.0 and 6.2. 

Our general conclusion is that the Wilson action at /3 = 6.0 has 0(a) errors of up to 20% 
compared to the continuum extrapolation. The non-perturbatively 0(a) improved theory 
still shows O(a^) effects of up to 10% at /5 = 6.0, except for the Ward identity quark masses 
where the effect is somewhat larger. If one wants to go to smaller values of /3, one probably 
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will have to reduce the O(a^) errors as well. Going to P = 6.2 reduces by a factor of 
almost two, bringing discretization errors down to 5% or less. To achieve a one percent 
accuracy would require calculations at several P values and an extrapolation to a = 0. 
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Tables 



/3 = 5.7 


csw = 1-0 


V 


K 




arup 


am AT 


aruao 


arua^ 




16^32 


0.1500 

0.1510 
0.1520 


0.5028(17) 

0.414(2) 

0.288(5) 


0.757(7) 
0.711(8) 
0.660(19) 


1.135(18) 
1.040(17) 
0.92(3) 


1.36(10) 
1.11(17) 


1.61(19) 
1.31(12) 
1.09(17) 


1.06(16) 
1.11(17) 
1.25(20) 


c. I. 


0.15280(14) 





0.605(24) 


0.797(49) 




0.70(36) 


1.31(30) 


^sw — 2.25 


V 






anip 


arriN 


arriao 






16^32 


0.1270 
0.1275 

0.1280 
0.1285 
0.1290 
0.1295 


0.841(3) 
0.791(4) 
0.736(3) 
0.672(5) 
0.607(7) 
0.519(11) 


1.087(8) 

1.053(10) 

1.022(11) 

0.988(9) 

0.955(8) 

0.922(16) 


1.588(12) 

1.518(23) 
1.453(18) 
1.399(24) 
1.320(20) 
1.23(3) 


1.64(10) 

1.56(7) 

1.50(7) 

1.57(14) 

1.59(14) 


1.56(7) 
1.51(5) 
1.46(4) 
1.41(6) 
1.34(8) 
1.28(10) 


1.48(13) 

1.53(9) 

1.42(7) 

1.39(7) 

1.33(11) 

1.33(16) 


c. I. 


0.13074(29) 





0.793(19) 


0.948(46) 


1.43(33) 


1.06(16) 


1.13(25) 



Table 1: Hadron masses at /5 = 5.7 for Sheikholeslami-Wohlert fermions with csw = 1 and 
2.25. In the bottom row we give Hc and the mass values extrapolated to the chiral limit. 
The numbers in roman {italic) are from three-parameter (two-parameter) fits. The errors 
are bootstrap errors. 
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/3 = 6.0 


Csw = 


V 






anip 


arriN 


arriao 


arria^ 


arrib^ 


16332 


0.1487 
0.1515 
0.1530 
0.1550 


0.6384(18) 
0.5037(8) 
0.4237(8) 
0.3009(10) 


0.683(2) 
0.5696(10) 
0.5080(11) 
0.4264(14) 


1.071(7) 
0.9019(17) 
0.7977(20) 
0.6517(30) 


0.885(19) 
0.817(7) 
0.763(11) 
0.735(15) 


0.933(13) 
0.851(7) 
0.797(6) 
0.717(12) 


0.940(19) 
0.849(13) 
0.809(7) 
0.736(9) 


24^32 


0.1550 
0.1558 
0.1563 


0.292(2) 
0.229(2) 
0.179(3) 


0.418(5) 
0.384(7) 
0.358(11) 


0.638(8) 

0.555(12) 

0.488(22) 


0.610(48) 
0.616(90) 
0.88(15) 


0.657(33) 
0.613(41) 
0.584(52) 


0.659(35) 
0.638(38) 
0.615(44) 


c. /. 


0.15713(3) 





0.327(6) 


0.412(16) 


0.658(19) 


0.632(14) 


0.650(13) 


Csw = 1-769 


V 






arrip 


arriN 


arriao 


arria^ 


arribj 


16^32 


0.1300 
0.1310 
0.1320 
0.1324 
0.1333 
0.1342 


0.707(2) 

0.627(2) 

0.545(5) 

0.5039(7) 

0.4122(8) 

0.2988(17) 


0.783(6) 

0.714(3) 

0.644(8) 

0.6157(16) 

0.5502(23) 

0.487(3) 


1.190(6) 

1.079(7) 

0.974(16) 

0.932(4) 

0.821(5) 

0.705(9) 


0.779(14) 
0.738(15) 
0.92(5) 


0.829(12) 

0.773(7) 

0.68(2) 


0.853(7) 

0.799(10) 

0.775(15) 


24^32 


0.1342 
0.1346 
0.1348 


0.3020(11) 
0.2388(14) 
0.194(4) 


0.491(3) 
0.467(6) 
0.448(13) 


0.686(7) 

0.626(10) 

0.593(19) 


0.82(3) 
1.00(8) 
1.52(20) 


0.715(19) 
0.684(26) 
0.664(34) 


0.758(16) 
0.745(20) 
0.736(29) 


c. /. 


0.13531(1) 





0.417(7) 


0.511(15) 


0.816(33) 


0.625(19) 


0.710(14) 



Table 2: Hadron masses aX (3 = 6.0 for Wilson fermions [csw = 0) and improved fermions 
{csw = 1-769). Otherwise the notation is the same as in table |I|. 
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P = 6.2 


Csw = 


V 






aiTip 


arriN 


arriao 


arua^ 




24^48 


0.1468 
0.1489 
0.1509 
0.1518 
0.1523 


0.5258(12) 
0.4148(13) 
0.2947(14) 
0.2299(15) 
0.1867(17) 


0.5585(16) 

0.4615(19) 

0.3672(27) 

0.326(4) 

0.307(6) 


0.872(5) 

0.720(6) 

0.560(10) 

0.487(12) 

0.448(14) 


0.685(8) 

0.589(8) 

0.507(14) 

0.474(20) 

0.479(30) 


0.700(21) 

0.624(9) 

0.536(13) 

0.509(16) 

0.492(17) 


0.695(21) 

0.626(9) 

0.540(13) 

0.519(17) 

0.511(21) 


c. /. 


0.15336(4) 





0.255(9) 


0.342(28) 


0.407(17) 


0.449(15) 


0.464(16) 


Csw = 1-614 


V 






arup 


arriN 


aniao 






24^48 


0.1321 
0.1333 
0.1344 
0.1349 
0.1352 


0.5179(7) 
0.4143(8) 
0.3046(9) 
0.2444(9) 
0.2016(11) 


0.5738(11) 
0.4850(15) 
0.4005(26) 
0.3626(43) 
0.3430(53) 


0.877(4) 

0.735(5) 

0.592(9) 

0.521(13) 

0.485(6) 


0.691(5) 

0.603(10) 

0.532(21) 

0.543(22) 

0.646(53) 


0.723(6) 

0.642(5) 

0.563(7) 

0.529(10) 

0.514(13) 


0.727(6) 

0.638(8) 

0.566(9) 

0.539(12) 

0.523(25) 


c. /. 


0.13589(2) 





0.287(9) 


0.378(18) 


0.460(21) 


0.460(9) 


0.465(12) 



Table 3: Hadron masses aX (3 = 6.2 for Wilson fermions {csw = 0) and improved fermions 
{csw = 1.614). Otherwise the notation is the same as in table |I|. 





Csw 


eq. (D 


eq. (0) 


eq. (D 








xVdof 










5.7 


1.0 


0.15305(5) 


3.0 


0.15274(15) 




0.15280(14) 




2.25 


0.13120(7) 


0.7 


0.13065(28) 


0.1 


0.13074(29) 


0.1 


6.0 





0.15695(1) 


17.5 


0.15726(5) 


8.6 


0.15713(3) 


6.6 


1.769 


0.13521(1) 


11.5 


0.13537(2) 


1.5 


0.13531(1) 


1.0 


6.2 





0.15308(1) 


30.8 


0.15361(8) 


0.7 


0.15336(4) 


0.0 


1.614 


0.13574(1) 


39.6 


0.13601(3) 


1.2 


0.13589(2) 


0.1 



Table 4: The critical values of Kc, of our data for the linear (eq. (P)), chiral (eq. (^) and 
phenomenological fit (eq. (||)) for the various csw parameters. 
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arUp 


arriN 


Lattice 


Reference 


6.30 


0.1400 


0.789(4) 


0.804(4) 




32^ X 48 




16 


1 


6.30 


0.1430 


0.646(6) 


0.670(5) 




32^ X 48 




m 


6.30 


0.1460 


0.4879(12) 


0.5188(18) 


0.8252(42) 


24^ X 32 




17 


1 


6.30 


0.1480 


0.382(4) 


0.429(4) 




32^ X 48 




16] 


6.30 


0.1485 


0.3480(14) 


0.3990(23) 


0.6340(47) 


24^ X 32 




17 


1 


6.30 


0.1498 


0.2631(19) 


0.3354(30) 


0.5215(67) 


24^ X 32 




17 


1 


6.30 


0.1500 


0.253(6) 


0.333(4) 




32^ X 48 




16| 


6.30 


0.1505 


0.2093(26) 


0.3012(40) 


0.4506(89) 


24^ X 32 




17 


1 


6.20 


0.1468 


0.5258(12) 


0.5585(16) 


0.872(5) 


24^ X 48 


this work 


6.20 


0.1489 


0.4148(13) 


0.4615(19) 


0.720(6) 


24^ X 48 


this work 


6.20 


0.1509 


0.2947(14) 


0.3672(27) 


0.560(10) 


24^ X 48 


this work 


6.20 


0.1510 


0.289(1) 


0.366(2) 


0.566(4) 


24^ X 64 




181 


6.20 


0.1515 


0.254(1) 


0.343(3) 


0.525(6) 


24^ X 64 




18 




6.20 


0.1518 


0.2299(15) 


0.326(4) 


0.487(12) 


24^ X 48 


this work 


6.20 


0.1520 


0.220(7) 


0.327(9) 


0.495(10) 


24^ X 48 




19| 


6.20 


0.1520 


0.215(1) 


0.321(5) 


0.48(1) 


24^ X 64 




18] 


6.20 


0.1523 


0.1867(17) 


0.307(6) 


0.448(14) 


24^ X 48 


this work 


6.20 


0.1526 


0.158(1) 


0.29(1) 


0.45(3) 


24^ X 64 




18] 


6.17 


0.1500 


0.3866(12) 


0.4458(18) 


0.6966(40) 


322 X 30 X 40 




2^ 


6.17 


0.1519 


0.2631(12) 


0.3572(26) 


0.5460(52) 


322 X 30 X 40 




201 


6.17 


0.1526 


0.2064(15) 


0.3245(39) 


0.4848(68) 


322 X 30 X 40 




2C| 


6.17 


0.1532 


0.1455(20) 


0.2965(88) 


0.4097(78) 


322 X 30 X 40 




2C 


1 



Table 5: World Wilson fermion masses above (3 = 6.0. 
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p 






arrip 


arriisi 


Lattice 


Reference 


D.U 


u 


145U 


U.8U6y(7j 


U.837U(9) 


1.3225(28) 


24^^ 


X 


54 


21 




o.U 


u 


1487 


U.6384(18j 


n POO f ^\ 

U. 683(2) 


l.U71(7) 


16'' 


X 


32 


this work 


o.U 


u 


1515 


U.5U37(8j 


U.5696(1U) 


U.9U19(17) 


16'' 


X 


32 


this work 


o.U 


u 


152U 


n A 'V70 { c\\ 

U.4772(9j 


n x:^ AOP ( A t:^\ 

U. 5486(15) 


f\ OP PC\ { AC\\ 

U. 8669(49) 


24^ 


X 


54 


|21 


1 


o.U 


u 


152U 


U. 474(1) 


U. 545(2) 


U. 861(5) 


18'' 


X 


32 


|17 


1 


D.U 


u 


153U 


U. 423(1) 


U.5U8(3) 


r\ on 1 { P\ 

U.8U1(6) 


18'' 


X 


64 


m 


1 


o.U 


u 


153U 


U. 4237(8) 


U.5U8U(11) 


U.7977(2U) 


16 


X 


32 


this work 


p. n 
O.U 


u 


153U 


U.4ZZ(i ) 


U.OUO(i ) 


U. ( oO(o ) 


32^ 


X 


64 


22 


1 


6.U 


u 


154U 


U.364(l) 


U.468(4) 


U.729(7) 


18^ 


X 


64 


18 




6.U 


u 


1545 


U.33U76(28) 


U.4425(1U) 


U. 6777(21) 


24^ 


X 


64 


|23 


1 


6.U 


u 


155U 


U.298(l) 


U.431(6) 


U.66(l) 


18^ 


X 


64 


|18 




6.U 


u 


155U 


U.3UU9(1U) 


U.4264(14) 


U.6517(3U) 


16^ 


X 


32 


this work 


6.U 


u 


155U 


U.29642(27) 


U.422U(12) 


U.6393(27) 


24^ 


X 


64 


|23 


1 


6.U 


u 


155U 


U.292(2) 


U.418(5) 


U.638(8) 


24^ 


X 


32 


this work 


6.U 


u 


155U 


U.2967(15) 


U.4218(42) 


U.644U(85) 


24^ 


X 


54 


|21 


1 


6.U 


u 


155U 


U.296(l) 


U.422(2) 


U.63U(5) 


32^ 


X 


64 


22 


1 


6.U 


u 


1555 


U.25864(33) 


U.4U16(17) 


U.6UU3(37) 


24^ 


X 


64 


23 


1 


6.U 


u 


1555 


U.2588(16) 


U.3982(61) 


U.6UU7(1U9) 


24^ 


X 


54 


21 




6.U 


u 


1558 


U.234(l) 


U.387(3) 


U.557(7) 


32^ 


X 


64 


22 


1 


6.U 


u 


1558 


U.229(2) 


U.384(7) 


U.555(12) 


24^ 


X 


32 


this work 


6.U 


u 


1563 


U. 1847(27) 


U.353(15) 


U.536(3U) 


24^ 


X 


54 


|21 


1 


6.U 


u 


1563 


U. 185(1) 


U.361(5) 


U.5U6(11) 


32^ 


X 


64 


22 




6.U 


u 


1563 


U. 179(3) 


U.358(ll) 


U.488(22) 


24^ 


X 


32 


this work 



Table 6: World Wilson fermion masses at (3 = 6.U. 
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/3 






arrip 


aniN 


Lattice 


Reference 


5.93 


0.1543 


0.4572(26) 


0.5527(40) 


0.8674(102) 


24^ X 36 


m 


5.93 


0.1560 


0.3573(19) 


0.4864(42) 


0.7448(99) 


24^ X 36 


m 


5.93 


0.1573 


0.2641(25) 


0.4369(48) 


0.6423(80) 


24^ X 36 


m 


5.93 


0.1581 


0.1885(31) 


0.4071(57) 


0.5652(92) 


24^ X 36 


m 


5.85 


0.1440 


1.0293(12) 


1.0598(15) 


1.6961(50) 


243 X 54 


m 


5.85 


0.1540 


0.6122(11) 


0.6931(27) 


1.1060(55) 


243 X 54 


m 


5.85 


0.1585 


0.3761(12) 


0.5294(69) 


0.815(13) 


243 X 54 


m 


5.85 


0.1585 


0.378(2) 


0.530(6) 


0.783(10) 


16^ X 32 


m 


5.85 


0.1595 


0.3088(14) 


0.4856(96) 


0.744(17) 


243 X 54 


m 


5.85 


0.1600 


0.2730(30) 


0.486(9) 


0.673(9) 


16^ X 32 


m 


5.85 


0.1605 


0.2226(21) 


0.434(20) 


0.683(48) 


243 X 54 


m 


5.70 


0.1600 


0.6905(31) 


0.8022(56) 


1.3124(135) 


243 X 32 


m 




5.70 


0.1600 


0.6873(24) 


0.8021(29) 


1.2900(60) 


16^ X 20 


m 




5.70 


0.1610 


0.6527(15) 


0.7842(26) 


1.263(5) 


12^ X 24 


m 




5.70 


0.1630 


0.5621(18) 


0.7232(35) 


1.153(6) 


12^ X 24 


M 


5.70 


0.1640 


0.5080(29) 


0.6822(38) 


1.0738(80) 


16^ X 20 


m 


5.70 


0.1650 


0.4604(22) 


0.6663(45) 


1.039(8) 


12^ X 24 


m 


5.70 


0.1650 


0.4589(22) 


0.6491(73) 


1.0301(104) 


243 X 32 


m 


5.70 


0.1663 


0.3829(26) 


0.6206(103) 


0.9421(131) 


243 X 32 


m 


5.70 


0.1665 


0.3674(39) 


0.6085(58) 


0.915(11) 


16^ X 20 


m 


5.70 


0.1670 


0.3302(30) 


0.6042(83) 


0.919(14) 


12^ X 24 




5.70 


0.1675 


0.2955(24) 


0.5912(125) 


0.8668(177) 


243 X 32 


^o| 





Table 7: World Wilson fermion masses below f3 = 6.0. 
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/5 




Reference 


ro/a 


Reference 




0.8 


0.0730(12) 


2§ 


16.7(4) 


Ih HII 


1.22(4) 


0.5 


0.1068(10) 


30| 








p. A 
0.4 


0.1215(12) 
U.iZio(Zo ) 


|29| 
1 Tm 


9.87(8) 
y. ( u(/4j 










Combined 


Q self's^ 


Combined 


i. iy ( (^ioj 


D.d 


n 1 on /I 1 1 \ 
U.ioy4(ii ) 


Interpolated 








6.2 


n 1 ^^1 n/'n\ 
U.iOiU(y ) 

n 1 Rn8('9'?^ 

U. iOUol^Zo ) 




r .00(4 j 


\'\^ \ 




0.1609(28) 


i33[ 










0.1610(8) 


Combined 


7.36(4) 


Combined 


1.185(9) 


6.17 


0.1677(8) 


Interpolated 








6.0 


0.2209(23) 
0.2154(50) 


m 


5.28(4) 
5.53(15) 




[0 




0.2182(21) 
0.2191(15) 


34 1 

Combined 


5.30(4) 


Combined 


1.161(12) 


5.93 


0.2536(29) 


Interpolated 








5.90 


0.2702(37) 


m 


4.62(11) 


[0 


1.25(3) 


5.85 


0.2986(27) 


Interpolated 








5.8 


0.3302(30) 


135[ 


3.63(5) 


m 


1.199(20) 


5.7 


0.4099(24) 


135[ 


2.86(5) 


m 


1.172(22) 


5.6 






2.29(6) 







5.5 






2.01(3) 








Table 8: The lattice spacing expressed in terms of the string tension K and the force param- 
eter tq. When several groups have computed these quantities, we have taken the weighted 
average, while we interpolate logarithmically whenever the values are not known. 
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/V,. 


6.40 


0.150759(145) 


6.30 


0.151774(36) 


6.20 


0.153374(17) 


6.17 


0.153838(37) 


6.00 


0.157211(8) 


5.93 


0.158985(73) 


5.85 


0.161716(23) 


5.70 


0.169313(72) 



Table 9: The critical values of Kc, for the Wilson world data. 



/3 = 6.0 





csw - 


= 


V 


K 


2 am 




0.1487 


0.2959(5) 


16^32 


0.1515 


0.1866(5) 


0.1530 


0.1321(5) 




0.1550 


0.0612(7) 





Csw — 


L.769 


V 


K 


2 am 




0.1300 


0.2836(3) 


16^32 


0.1310 
0.1324 


0.2279(3) 
0.15231(10) 




0.1333 


0.10380(11) 




0.1342 


0.0553(2) 




0.1342 


0.0551(3) 


24^32 


0.1346 


0.0330(3) 




0.1348 


0.0214(4) 



Table 10: The bare quark masses m for Wilson fermions {csw — 0) and improved fermions 
{csw — 1.769) at P — 6.0. The errors are bootstrap errors. 
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csw - 


-- 


V 


K 


2 am 




0.1468 


0.2474(3) 




0.1489 


0.1616(3) 


24^48 


0.1509 


0.0845(3) 




0.1518 


0.0514(3) 




0.1523 


0.0336(3) 



csw = 1-614 


V 


K 


2a7fi 




0.1321 


0.2161887) 




0.1333 


0.14585(7) 


24^48 


0.1344 


0.08185(7) 




0.1349 


0.05283(8) 




0.1352 


0.03538(9) 



Table 11: The bare quark masses m for Wilson fermions {csw — 0) and improved fermions 
{csw = 1.614) at /3 = 6.0. The errors are bootstrap errors. 





Csw 




m 


MS 
s 


Ward 


Standard 


Ward 


Standard 


6.0 





4.40 ±0.17 


6.47 ±0.20 


105.0 ±4.5 


141.8 ±6.0 


6.2 





4.73 ±0.14 


6.39 ±0.25 


108.5 ±4.2 


138.8 ±7.4 


6.0 


1.769 


4.02 ±0.10 


4.94 ±0.93 


92.8 ±2.9 


109.4 ±3.9 


6.2 


1.614 


4.47 ±0.06 


5.09 ±0.16 


101.6 ± 1.7 


111.7±4.7 


oo 




5.00 ±0.18 


5.27 ±0.36 


111.9 ±5.0 


114.4 ± 11.1 



Table 12: Our results of the renormalized quark masses m^'^(2GeV) in MeV for improved 
and Wilson fermions, together with the extrapolation to the continuum limit {(3 — oo). 
The continuum numbers refer to improved fermions. We give the results for both the Ward 
identity and the standard method. 
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/3 = 6.0 


csw = 


V 


K 










f(0) 
J ai 


16^32 


0.1487 
0.1515 
0.1530 

0.1550 


0.136(2) 
0.122(2) 
0.113(2) 

0.098(2) 




0.305(5) 
0.364(5) 
0.397(7) 

0.159(9) 




0.161(5) 
0.207(3) 
0.231(3) 

0.262(4) 


Csw = 1-769 


V 


/V 










fm 


16^32 


0.1300 
0.1310 
0.1324 
0.1333 

0.1312 


0.1341(15) 
0.1295(15) 
0.1204(8) 
0.1128(9) 

0.1037(8) 


1.792(7) 
1.698(8) 
1.599(5) 
1.541(3) 
1.511(6) 


0.209(9) 
0.228(3) 
0.261(2) 
0.288(3) 

0.323(3) 


0.670(2) 
0.588(2) 
0.4823(12) 
0.4147(15) 

0.353(2) 


0.131(7) 
0.153(13) 
0.172(7) 
0.202(16) 

0.208(16) 


24^32 


0.1342 
0.1346 
0.1348 


0.105(2) 
0.101(2) 
0.100(3) 


1.521(16) 

1.58(3) 

1.62(5) 


0.330(7) 
0.352(7) 
0.352(15) 


0.348(3) 
0.327(6) 
0.324(19) 


0.212(10) 
0.225(14) 
0.25(2) 



Table 13: The various contributions to the decay constants fp and fa^ ai (3 — 6.0. 
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p = 6.2 



Csw = 


V 












m 

J ai 


24^48 


0.1468 
0.1489 
0.1509 
0.1518 
0.1523 


0.1025(19) 
0.0930(17) 
0.0798(14) 
0.0719(13) 
0.0669(14) 




0.268(5) 
0.315(6) 
0.376(7) 
0.412(9) 
0.438(10) 




0.127(9) 
0.180(4) 
0.230(4) 
0.261(4) 
0.276(5) 


csw = 1.614 


V 












m 

J ai 


24^48 


0.1321 
0.1333 
0.1344 
0.1349 
0.1352 


0.0985(11) 

0.0913(11) 

0.0818(10) 

0.0758(9) 

0.072(3) 


1.297(3) 
1.198(4) 
1.133(2) 
1.118(7) 
1.131(11) 


0.211(3) 
0.243(3) 
0.283(4) 
0.308(5) 
0.327(16) 


0.4637(8) 

0.3719(10) 

0.2900(15) 

0.255(2) 

0.235(4) 


0.133(2) 
0.167(2) 
0.204(3) 
0.226(3) 
0.241(11) 



Table 14: The same as table 13 but for P = 6.2. 





Csw 




fai 








6.0 
6.2 






0.0569(77) 
0.0423(36) 


0.2240(76) 
0.2429(48) 


0.295(11) 
0.2971(71) 


0.0732(26) 
0.0537(11) 


0.2385(24) 
0.2345(28) 


6.0 
6.2 


1.769 
1.614 


0.0627(20) 
0.0462(38) 


0.195(10) 
0.2130(45) 


0.2664(39) 
0.2726(74) 


0.0721(8) 
0.0557(14) 


0.2020(13) 
0.2149(19) 



Table 15: The decay constants fa^, fp extrapolated to the chiral limit, as well as 
taken at the physical quark mass. 
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Figures 



0.15 



0.14 - 



0.13 - 



0.12 




Figure 1: The critical values of n for improved Wilson fermions as a function of from this 
work (o) and ref. (x). The curve is the tadpole improved result given in eq. ([T0|). 
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0.2 0.4 0.6 0.2 0.4 0.6 



(amj2 (am^)2 

Figure 2: Fits and chiral extrapolations of p and nucleon masses for improved (•) and Wilson 
fermions (o). 
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Figure 3: APE plot a,t P — 6.0 for improved (•) and Wilson fermions (o) compared with 
the physical mass ratio (^) at the physical quark mass and in the heavy quark limit. The 
solid lines are from the mass fits described in the text. 



33 



1.6 




0.2 0.4 0.6 0.8 1 

(mym^)2 



Figure 4: The same as fig. but for (3 = 6.2. Tlie liatclied bars indicate tlie errors in 
cliiral limit. 
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m/VK 




1 I — ' — ' — ' — ' — ^ — ' — ' — ' — ' — ^ — ' — ' — ' — ' — ^ — ' — ' — ' — ' — ^ — ' — ' — ' — ' — I 
0.1 0.2 0.3 0.4 0.5 

a Vk 



Figure 5: The ratio nip/ \/K as a function of the lattice spacing for improved (•) and Wilson 
fermions (o). The solid lines are from a simultaneous linear plus quadratic fit to the Wilson 
data and a quadratic fit to the improved data. This is compared with the experimental value 
(*) using y/K = 427 MeV. 



35 




5 10 15 



mVK 

Figure 6: The vector-pseudoscalar mass splitting as a function of the quark mass. Open 

symbols correspond to (3 = 6.0, solid symbols to f3 = 6.2. This is compared with the 
physical p-vr (^) and K*-K (octogon) mass splitting. The curves are from the mass fits. 
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a Vk 
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0.5 



Figure 7: The critical value of /t as a function of the lattice spacing for Wilson fermions 
(□). The dashed curve is the prediction of tadpole improved perturbation theory. This is 
compared with the results for improved fermions from fig. ||. 
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Figure 8: The bare mass m from the Ward identity method at f3 = 6.0 for Wilson (□) and 
improved fermions (■) on the 16'^32 lattice. The errors are bootstrap errors. 
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Figure 9: The same as fig. H but for (3 = 6.2 for Wilson (□) and improved fermions 
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Figure 10: The light quark mass as a function of the lattice spacing for improved 

fermions using the Ward identity (•) and standard method (■). This is compared with the 
Wilson result for the Ward identity (o) and standard method (□). The curves are quadratic 
extrapolations to the continuum limit. 
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Figure 12: The light quark mass m^^^ for improved fermions from fig. [T0| compared with the 
world Wilson masses (A) compiled in 
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Figure 13: The strange quark mass m^^^ for improved fermions from fig. |TT] compared with 
the world Wilson masses (A) compiled in These authors use the 0(1020) meson to 

determine the strange quark mass. 
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Figure 14: Fits and chiral extrapolations of the decay constants and fp for improved (•) 
and Wilson (o) fermions. 
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Figure 15: The decay constants and fp as a function of the lattice spacing for improved 
(•) and Wilson (o) fermions together with the experimental values (^). The errors on fp 
for improved fermions are statistical only. Our results for are compared with the Wilson 
results of ref. 1151 



45 




Figure 16: The decay constants fx and fx* as a function of the lattice spacing for improved 
(•) and Wilson (o) fermions. The errors on fx* for improved fermions are statistical only. 

The solid lines in the fx* figure are from a simultaneous linear fit to the Wilson data and a 
quadratic fit to the improved data. Our results for are compared with the experimental 
value (^). 
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